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ABSTRACT 


A  summary  of  the  Weibull  distribution  and  three  repre 
sentative  sampling  prooedures  to  determine  point  and  con¬ 
fidence  interval  estimates  of  the  parameters  that  ooour 
in  the  functional  form  of  the  Weibull  distribution  are 
presented.  Following  this,  a  model,  assuming  the  Weibull 
distribution,  is  proposed  which  could  have  possible  ap¬ 
plications  in  analyzing  Polaris  Missile  System  Trouble 
Failure  Reports  to  determine  point  estimates  of  the 
reliability  of  the  Polaris  Missile  components, 

I  am  indebted  to  Professor  P,  W,  Zehna  for  his  con¬ 
tinued  patience,  encouragement,  and  most  capable  guidance 
while  acting  as  faculty  advisor,  I  also  wish  to  thank 
Professor  W,  M.  Woods  for  his  valuable  assistance  as 
second  reader. 
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CHAPTER  I 


INTRODUCTION 


This  thesis  is  concerned  with  summarizing  the  Weibull 
distribution  and  providing  several  possible  methods  of 
analyzing  Polaris  Missile  System  Trouble  Failure  Reports 
to  determine  the  reliability  of  Polaris  Missile  com¬ 
ponents  when  the  Weibull  distribution  is  assumed  to 
characterize  the  failure  times  of  the  components* 

The  Weibull  distribution,  which  was  first  proposed 
by  Waloddi  Weibull  M  in  1951#  is  a  member  of  a  class  of 
distributions  which  characterize  wear  out  failure;  the 
longer  an  item  has  been  used  the  greater  the  probability 
of  failure*  Wear  out  failure  is  associated  with  gradual 
depletion  of  material,  fatigue,  or  accumulated  shocks 
and  other  factors* 

The  key  assumption  underlying  many  of  the  statistical 
procedures  in  current  use  for  evaluating  the  reliability 
of  components  or  equipment  is  that  the  failure  times 

follow  the  exponential  distribution.  However,  in  many 

♦ 

practical  situations  sufficient  data  is  not  usually 
available  to  Justify  the  assumption  that  the  exponential 
failure  law  truly  characterizes  the  time  to  failure 
distribution  of  the  equipment  or  components  being  tested. 
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The  Weibull  distribution  which  includes  the  exponential 
distribution  as  a  special  case  is  a  more  general  assump¬ 
tion  about  the  time  to  failure  law  of  components  or 
equipment*  Kao  £2^}  has  shown  that  the  failure  age 
distribution  of  electron  tubes  is  better  chacterized 
by  the  Weibull  distribution  with  shape  parameter  ap¬ 
proximately  equal  to  1.7  than  by  the  exponential  dis¬ 
tribution.  Also  when  the  exponential  distribution  is 
assumed  to  be  the  time  to  failure  law,  M.  Zelen  and 
M.  C.  Dannemiller  [3]  demonstrate  that  four  commonly 
employed  sampling  procedures  are  not  robust*  with  re¬ 
spect  to  Weibull  alternatives. 

Chapter  II  summarizes  the  distribution  theory 
associated  with  the  Weibull  distribution  and  defines 
terms  to  be  used  in  the  remainder  of  this  thesis. 

Chapter  III  outlines  sampling  plans  which  can  be  em¬ 
ployed  to  obtain  estimates  of  the  parameters  from  ex¬ 
perimental'  measurements.  The  various  estimates  to 

follow  are  derived  under  each  of  these  plans.  Chapter 

\ 

IV  proposes  three  possible  methods  of  using  Polaris 
Missile  Trouble  Failure  Reports  for  determining  the 
reliability  of  Polaris  system  components. 

*  Sampling  procedures  which  are  not  very  sensitive  to 
departures  from  basic  assumptions  are  called  robust. 


2 


CHAPTER  II 


SUMMARY  OP  THE  WEIBULL  DISTRIBUTION 


Let  X  be  a  random  variable  which  may  take  on  any 
value  in  the  oontinuous  scale  from  m  to  infinity  where 
m  is  a  non-negative  number.  The  three  parameter  Weibull 
cumulative  distribution  function  (b.d.f,)  is  defined  as 


(1).  F(x) 


r  -(x-m)** 
a 

-  1-e  for  x  t  m,  where 

0  x  <  m 


a  -  scale  parameter,  a  >  0 
b  -  shape  parameter,  b  >  0 
m  x  location  parameter*,  m  i  0 
The  probability  density  function  is  defined  as  the 
first  derivative  of  F(x)  with  respeot  to  x,  and  is  given  by. 


(2).  '  dP(x)  =  f(x)  _  < 

dx 


b(x-m) 


b-1 


-(x-m) 

a 


x  2s  m 


a 


V.0 


x  <  m 


Weibull  fi 1  demonstrated  that  yield  strength  of 
Bofors  steel,  size  distribution  of  fly  ash,  and  fiber 


*  xQ,  m,  and  Xm,  respectively  were  the  notations  used 
by  Weibull  QT]  • 
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strength  of  Indian  cotton  are  distributed  in  good  agree¬ 
ment  with  the  cumulative  distribution  function,  equation 


(1).  Kao  M  has  shown  that  the  cumulative  distribution 


function,  equation  (1),  qualifies  as  a  failure  distribu¬ 
tion  and  has  successfully  employed  the  Weibull  distribu¬ 
tion  with  location  parameter,  m,  equal  to  zero  to  char¬ 
acterize  the  time  to  failure  of  electron  tubes*  For  the 

purposes  of  this  thesis  and  the  sampling  plans  which 

* 

follow  in  Chapter  III,  the  location  parameter,  m,  will 
be  assumed  to  be  zero*  This  is  a  reasonable  assumption 
since  if  an  item  is  placed  in  operation  at  time  zero,  it 
is  exposed  to  the  risk  of  failure  from  the  time  it  is 
first  put  into  use. 

For  completness  of  presentation,  the  moments  and 
characteristic  function  of  the  three  parameter  Weibull 
distribution  are  listed  below,  but  a  detailed  discussion 
will  only  be  presented  for  the  two  parameter  sub-family 
having  msO.  Given  f(x)  as  the  above  probability  den¬ 
sity  function,  the  nth  moment  is  defined  by  the  formula 


-(x-m)k 
b-1  a 


E(Xn)  =  )  xn  b(x-m) 


e 


a 

and  is  given  by 


n 


k 


k 


4 

The  wean  and  variance,  denoted  by  X  and  (T  respectively, 
are  easily  determined  from  the  above  relation  to  be 

1 

ET  - .  .  l 

X  s  m+  a  J  1  (b  +  1) 

and 

(t  +  1)  (E  +  1)J 

respectively.  The  function  T<  • )  is  the  gamma 
function  which  is  defined  for  every  v  ?0  by 

_  r-°o 

p(v)  =  }  x*”1  e“x  dx 

Suppose  now  that  T  is  a  random  variable  which  may 
take  on  any  value  in  the  continuous  time  scale  frota  zero 
to  infinity,  and  the  random  variable  T  represents  the 
total  time  to  failure  of  an  item.  Let  us  further  assume 
that  the  two  parameter  Weibull  probability  density  func¬ 
tion,  f(t),  specifies  the  probability  law  of  the  random 
variable  T  where 

(3)«  f(t)  =  (btb**Va)  exp  (-tb/a),  t  *5.  0;  a,b,  >0 
—  0  elsewhere 

Accordingly,  if  an  item  is  placed  in  operation  at  time 
zero,  the  probability  that  the  total  life  of  the  item 
is  less  than  or  equal  to  t  is  expressed  as 
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P(T  5  t)  = 


Therefor©  the  cumulative  distribution  function  of  the 
random  variable  T  is  given  by 

(If.)*  F(t)  =  1  -  exp  («tb/a) 

The  reliability  function,  denoted  by  R(t),  of 
an  item  is  also  of  particular  interest  since  it  is 
desirable  to  know  the  probability  that  an  item  or 
component  does  not  fail  in  the  time  interval  (0,t). 

Thus  the  reliability  function  is  the  probability 
that  the  time  to  failure  of  an  item  is  at  least  t 
time"  units  and  is  expressed  by 

R(t )  s  P(T  >  t)  =  (btk“Va)exP(-tb/£L)dt,  for  t>0 

which  upon  integration  becomes 

% 

\ 

($).  R(t )  3  exp  (-tb/a)  i  t  ?0 
Note  that  the  reliability  function  is  merely  the  unity 
complement  of  the  cumulative  distribution  function, 
equation  (If.)* 
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The  nth  moment  of  the  random  variable  T  is  defined 
by  the  formula 

,  E(Tn)  =  ^tn(btb“Va)  exp  ( -tb/a )  dt 
and  is  given  by 

n 

(6).  E(Tn)  =  a  ]""*  (K  1) 

The  characteristic  function  of  the  random  var¬ 
iable  T  is 

$(u)  s  E(e*uT)  5  \  eiut  b  tb~1  e“(tbA)  dt 

Jo  a 

which  upon  integration  becomes 

(j>(u)  -  T1  (£+>1) 

The  mean  time  to  failure  of  the  random  variable  T, 
denoted  by  X  ,  is  defined  as  the  first  moment  about 
zero  and  Is  easily  determined  from  equation  (6)  to  be 

1 

F  *T  1 

(7).  |  (E+l) 

In  a  similar  manner,  the  variance  denoted  by  O'* of  the 
random  variable  T  is  determined  by  the  formula. 


t 
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E(T2)  -  E(T)2,  to  be 

<rA=  *  (i  *•'  i)  -~T 

Another  function  of  a  time  to  failure  probability 
distribution  which  provides  a  great  deal  of  insight 
into  the  probability  law  is  the  hazard  function  or 
instantaneous  failure  rate  which  will  be  denoted  by 
Z(t).  The  instantaneous  failure  rate  or  hazard 
function  is  defined  by 

where  f(t)  and  R(t)  are  the  density  and  reliability 
functions  respectively.  Substitution  of  f(t)  and 
R(t)  of  equations  (3)  and  (5)  in  the  above  express¬ 
ion  , we  have 

(8).  Z(t)  —  bt^"Va  for  tsrO 

It  is  interesting  to  note  that,  for  values  of  the  shape 
parameter  greater  than  one,  the  instantaneous  failure 
rate  is  an  increasing  function  of  time,  therefore  indi¬ 
cating  that  wear  out  occurs  during  the  life  of  the  item. 
But  in  the  special  case  where  the  shape  parameter  equals 
ofae  and  the  Weibull  distribution  simplifies  to  the  ex¬ 
ponential  distribution,  it  is  observed  that  the  instan- 
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taneous  failure  rate  is  a  oonstant  equal  to  1/a,  inde¬ 
pendent  of  age. 

The  effect  of  varying  the  shape  and  scale  para¬ 
meters  can  most  readily  be  illustrated  by  graphs  of 
the  probability  density  function  because  the  prob¬ 
ability  that  the  random  variable  T  will  lie  in  the 
interval  (t,  t-*-h)  can  be  interpreted  geometrically 
as  the  area  under  the  probability  density  function  in 
the  interval  from  t  to  t  +  h.  Figure  one  illustrates 
the  probability  density  function  for  several  values  of 
the  shape  parameter,  b,  and  with  the  scale  parameter, 
a,  equal  to  one.  The  effect  of  changing  the  value  of 
the  scale  parameter  is  to  merely  squeeze  or  broaden 
the  graph  of  the  probability  density  function,  i.e,, 
an  increase  in  the  value  of  the  scale  parameter 
would  broaden  while  a  decrease  in  the  value  of  the 
scale  parameter  would  squeeze  the  graph  of  the  prob¬ 
ability  density  function  together.  The  first  deriva¬ 
tive  of  f(t)  with  respect  to  t  is 

-tb 

f»(t)  -  b  tb“2  e  j~-b  tb  +  (b-l)J 

which  shows  that  f(t)  has  a  maximum  at 
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if  b  is  greater  than  one.  But  if  the  shape  parameter 
has  a  value  greater  than  zero  but  less  than  or  equal 
to  one,  the  probability  density  function  is  a  mono- 
tonically  decreasing  function  of  t.  Thus  increas¬ 
ing  values  of  the  shape  parameter  squeezes  the  pro¬ 
bability  density  function  together*  This  can  be 
interpreted  intuitively  to  mean  that  the  faster  an 
item  wears  out,  the  greater  the  value  of  the  shape 
parameter,  b*  Figure  2  illustrates  the  reliability 
function,  R(t),  for  several  values  of  the  shape 
parameter  and  with  the  scale  parameter,  a,  equal 
to  one.  Notice  for  fixed  values  of  the  scale  para¬ 
meter,  all  the  curves  intersect  at  t  equal  to  one, 
regardless  of  the  value  of  the  shape  parameter. 

Thus  far,  consideration  has  only  been  given  to  a 
single  item.  It  would  be  desirable  to  be  able  to  con¬ 
sider  a  component  or  complex  item  consisting  of  M 
individuals  or  single  items.  Instead  of  assuming  that 
the  items  are  independent,  it  will  be  assumed  that  the 
M  items  of  the  complex  item  are  quasi- independent  in 
the  following  sense.  If  the  life  of  the  ith  item  is 
a  random  variable  T^  with  corresponding  density  func¬ 
tion  F^(t)  then  let  T  be  a  random  variable  whose 
cumulative  distribution  function,  G(t),  is  determined 
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R(t> 


>  / 
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Figure  2.  Graph  of  reliability  function  for  several  values  of  shape  parameter,  b,  and 

scale  parameter,  a,  equal  to  one. 


by  the  equation 


l-G(t)  =  P(TX  >  t)P(Ta^  t)  .  .  .  P(TM>t) 


for  t  greater  than  or  equal  to  zero 
the  above  equation  for  G(t)  yields 


Solving 


T  i1 

ir  a: 


t*i 

a. 


(9).  G(t)  =•  1-e  ,  t^O 

The  first  derivative  of  G(t)  with  respect  to  t  de¬ 
noted  by  g(t)  gives 


(10). 


dG(t)  -  g(t)s  V  bit 

dt  4r~ I 


b;-l 


Equations  (9)  and  (10)  above  will  define  the  pro¬ 
bability  density  and  cumulative  distribution  func¬ 
tion  respectively  of  a  complex  ,  i.e.,  a 

component  consisting  of  more  than  one  item  which  are 

quasi-independent.  The  reliability  function  of  a 

\ 

complex  item  is  determined  by  the  probability  state¬ 
ment,  P(T>t),  to  be 


-r 


A|  * 

?! 

7T/  ai 


(11).  R(t)  —  \  g(t)  dtse  A  ,  t>0 

ahd  the  hazard  function  is  determined  from  the  re¬ 


lation  Z  ( t  )  ««  g(t) 
l- 


to  be 


GTtT 


) 
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(12).  2(t)o^b1tbi"1/»1 

The  formula  for  computing  the  nth  moment  of  a  com¬ 
plex  item  is  then  given  by 

•  oo 


Ed11)  =■  ^  tn  'y'  (bjtkj^/a^)  exp  (- ^t^/a^Jdt 


M 

r 

'O  l,t  - 

which  can  be  integrated  by  parts  to  yield 


5-®  M 

,  ntn-1  exp  (-  ^  tbi/a. )dt 

«  £7 

Again  the  mean  and  variance  will  be  denoted  by  A  and  (T 
respectively  and  are  determined  by  equation  (13)  to  be 


(111.) 


r°°  m 

•  —  j  exp  (-  2  tbi/ai)<lt 

*4 o  x?7  i 


A  C  **  V\  a. 

(15>K  (T  *  )  2t  exp  (-  5Z  t®1/61!)^  "  A 

Equations  t( Xip)  and  (l£)  ©an  be  evaluated  by  numerical 
methods. 

If  the  Weibull  family  of  probability  distributions 
is  assumed  to  characterize  the  time  to  failure  of  an 
item  or  component  consisting  of  more  than  one  item 
where  the  various  members  of  the  family  differ  only 
by  the  values  of  the  parameters  occuring  in  the  func¬ 
tional  form  of  the  distribution,  it  is  dear  that  the 
reliability,  hazard,  density,  and  cumulative  distri- 


but ion  functions,  along  with  the  various  moments  are 
functions  of  the  parameters,  a  and  b.  Consequently, 
useful  probability  statements  about  these  items  can 
only  be  made  to  the  extent  that  estimates  of  these 
parameters  may  be  given.  At  the  present  time,  point 
and  confidence  interval  estimation  of  parameters  are 
most  frequently  employed. 

The  maximum  likelihood  estimates  of  the  para¬ 
meters,  a  and  b,  will  be  derived  for  each  of  the 
sampling  procedures  presented  in  chapters  III  and 
IV,  although  confidence  interval  estimates  will  only 
be  enumerated  for  the  scale  parameter,  a,  when  the 
shape  parameter,  b,  is  assumed  to  be  known  in  one  of 
the  sampling  plans.  A  thorough  search  of  recent  tech¬ 
nical  publications  indicates  that  confidence  interval 
estimates,  when  both  parameters  are  unknown,  have  not 
been  derived  and  the  derivation  of  these  confidence 
interval  estimates  is  beyond  the  scope  of  this  thesis. 
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CHAPTER  III 
SAMPLING  PLANS 

In  this  chapter,  several  methods  of  obtaining 
experimental  data  to  determine  the  maximum  likelihood 
estimates  (M.L.E.)  of  the  parameters  that  appear  in  the 
functional  form  of  the  two  parameter  Weibull  probability 
distribution  will  be  discussed.  Experimental  data  can 
be  obtained  in  a  variety  of  ways,  and  each  particular 
method  of  obtaining  experimental  measurements  will  be 
defined  as  a  sampling  plan.  For  each  of  the  sampling 
plans  presented,  it  is  assumed  that  the  random  var¬ 
iable  T  characterizes  the  time  to  failure  of  an  item 
and  the  density  function,  f(t),  defined  by  equation 
(3),  specifies  the  probability  law  of  T.  It  is  fur¬ 
ther  assumed  that  it  is  possible  to  make  repeated 
measurements,  denoted  by  T^,  Tg#  ♦  •  •  ,  Tn,  of  the 
random  variable  T,  and  the  measurements,  T^,  •  .  •  ,Tn, 
are  independent  identically  distributed  random  var¬ 
iables,  each  with  the  same  distribution  as  the  ran¬ 
dom  variable  T.  Accordingly,  after  selecting  a  ran¬ 
dom  sample  and  performing  an  experiment,  the  observed 
values  of  the  random  variable  T  are  designated  by 
t^,  tg,  ...  »t^« 


) 
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In  the  case  of  a  complex  item  composed  of  M 
individual  items,  a  random  sample  of  size  n^  will 
be  taken  from  the  item  of  the  ith  type,  for 
i  s  1,  2,  •  •  •  ,M,  Again  the  random  variables 
Tij'  for  f 1x6(1  1  and  J  =  2,  ,  .  .  ,  ni,  are 

independent  identically  distributed  random  var¬ 
iables*  Also  after  performing  ah  experiment,  the 
observed  values  of  the  random  variable  T^  will  be 
denoted  by  t^,  tj^,  •  *  t^n^.  Each  of  the 

sampling  plans  presented  are  based  upon  random 
samples  as  defined  in  the  above  paragraph  and  are 
restatements  of  the  sampling  plans  contained  in 
reference  . 

The  method  of  maximum  likelihood  is  utilized 
to  obtain  estimates  of  the  parameters,  because  this 
procedure  generally  tends  to  have  desirable  proper¬ 
ties  not  always  shared  by  other  methods.  In  brief, 
the  method  of  maximum  likelihood  estimation  consists 
of  determining  the  values  of  the  parameters  that  max¬ 
imize  the  likelihood  function,  denoted  by  L.  The 
likelihood  function,  L,  is  a  function  of  the  para¬ 
meters  in  the  joint  density  with  fixed  observations. 


) 
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This  maximizing  process  is  generally  accompolished  by 
determining  the  values  of  the  parameters  that  maximize 
the  natural  logarithm  of  L,  since  In  L  is  an  increas¬ 
ing  function  of  L.  The  functional  form  of  L  will  vary 
in  general  with  each  sampling  plan.  Lloyd  and  Lipow 


and  Mood  [7]  provide  a  detailed  discussion  of 


the  maximum  likelihood  method. 

The  maximum  likelihood  estimates  of  the  para- 

A  A 

meters,  a  and.  b,  will  be  denoted  by  a  and  b  res¬ 
pectively  in  each  of  the  sampling  plans  discussed 
in  this  paragraph.  These  sampling  plans  are  out¬ 
lined  as  follows: 

SAMPLING  PLAN  I  -  Fixed  sample  size. 

Single  item  case  -  A  random  sample  of  size  n  is  selec¬ 
ted  and  each  member  of  the  sample  is  tested  until  it 
fails.  The  total  time  to  failure  of  each  member  of 
the  sample  is  observed  and  recorded.  If  the  obser¬ 
ved  failure  times  of  the  random  variable  T  are 
t^,  tg»  ...»  tn»  then  the  likelihood  function 
is  given  by 


L  5=  "77"  f(ti) 

J-/ 


Upon  substituting  equation  (3)  for  f(t^)  in  the 
above  relation  and  taking  the  natural  logarithm 


18  . 


of  L,  we  have 


n  n 

In  L:n  In  b  •  n  In  a  +  ^  In  -  ^t,b/a 

A*l  1 

Therefore,  if  both  the  parameters,  a  and  b,  are  un¬ 
known,  the  M.L.E.  of  these  parameters  are  given  by 

n  S' 

(16) .  a*  s  t^/n 

^  ^  A  yy 

(17) .  b  ^  na/(  %  tib  la  t4  -  a  2lntJ 

I7»  x 

respectively.  Lloyd  and  Lipow  ^6^J  outline  a  pro¬ 
cess  of  iteration  which  will  solve  equations  (16) 

,  .  A  a 

and  (17)  for  a  and  b.  In  the  event  that  the  para¬ 
meter,  b,  is  assumed  to  be  known,  the  M.L.E.  of  the 
parameter,  a,  under  this  sampling  plan  is  given  by 

(18) .  a“  a  2  t4b/n 

In  this  situation,  it  is  also  observed  that  a  con¬ 
fidence  interval  estimate  for  the , scale  parameter 
can  easily  be  determined  using  the  statistic, 

'  (2  *£  t±b)/a 

which  has  a  Chi-square  distribution  with  2n  degrees 
of  freedom.  Complete  details  of  confidence  inter¬ 
val  estimates  can  be  found  in  Mood  M  • 
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Complex  Item  case  -  Assuming  that  the  complex  item  is 
composed  of  M  individual  items,  the  procedure  is  to 
select  a  random  sample  of  size  n^  from  the  item  of 
the  ith  type,  for  is  1,  2,  *  .  .  ,  M.  Each  member 
of  the  ith  random  sample  is  tested  until  it  fails, 
and  the  total  time  to  failure  is  observed  and  re¬ 
corded  for  each  member  of  the  ith  sample.  If  the 
observed  failure  times  of  the  random  variable 
are  t^,  t^,.  •  •  •  *  t^n^,  then  the  likelihood 
function  for  the  item  of  the  ith  type  is  given  by 


m 


Li  = 


for  i:l,  2,  .  •  •  ,M.  Therefore,  upon  subs¬ 
tituting  equation  (3)  for  f(t^j)  in  the  above  re¬ 
lation,  we  obtain 


L*  as  ~7f(*>i  exp  (-t^i/a*) 

Accordingly,  the  M.L.E.  of  the  scale  and  shape  para- 
meters  for  the  item  of  the  ith  type  are  determined, 
in  the  same  manner  as  for  a  single  item,  to  be 

(19).  -  Z?  t1Jbl/n1 

and 

:  <2°)*  ln 

-ai  fr ln  V 
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respectively,  for  i  =  1,  2,  .  .  .  ,  M.  Again,  a 
process  of  iteration  oan  be  employed  to  solve  for 


the  mamimum  likelihood  estimates.  In  the  situation 
where  all  of  the  shape  parameters  are  assumed  to  be 
known,  the  M.L.E.  of  a^  under  this  sampling  plan  i3 
given  by 


for  i  s  1,  2,  •  ,  ,  ,  M* 

SAMPLING  PLAN  II  -  Fixed  sample  size  with  item 
truncation. 

Single  item  case  -  A  random  sample  of  size  n  i3 
placed  on  te^t  simultaneously,  and  the  life  test 
is  terminated  when  r  members  of  the  sample  fail. 
The  number  of  failures  r  is  fixed  in  advance  of 
the  experiment,  and  r  is  an  integer  greater  than 
zero  but  less  than  or  equal  to  n.  Using  thi3 


sampling  plan,  the  first  r  failures  times  of  the 
random  variable  T  are  observed  and  recorded  as 
the  failures  occur.  Then,  if  ,  •  .  , 


i  tr<  oo  are  the  first  r  failures,  Halperin  QO 


indicates  that  the  likelihood  function  is  given  by 
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where  L  is  merely  the  density  funotion  of  the  first 
r  order  statistics.  After  substituting  F(t)  and  f(t) 
of  equation  (3)  and  (4)  in  the  above  equation,  the 
M.L.S.  of  the  parameters,  a  and  b,  are  determined  in 
the  usual  manner  to  be  given  by 


Again,  the  values  of  a  and  b  which  solve  equations 
(22)  and  (23)  can  be  determined  by  an  iteration 
process.  In  the  event  that  the  shape  parameter 
is  assumed  to  be  known,  the  M.L.E.  of  the  para¬ 
meter,  a,  under  this  sampling  plan  is  given  by 

(24). ,»  a  -  (  2?  tib  *  (n-r)trb)/r 

Complex  item  case  -  Assuming  that  the  complex  item 
is  composed  of  M  individual  items,  the  procedure 
is  to  select  a  random  sample  of  size  n,,  from  the 
item  of  the  ith  type,  for  i  -  1,  2,  ...»  M.  The 
n^  members  of  the  ith  random  sample  are  placed  on  test 


1 
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simultaneously  and  the  ith  life  test  is  terminated 
when  members  of  the  ith  sample  fall.  As  before, 
the  number  of  failures  ri  i3  fixed  in  advance  of 
each  experiment,  and  r^  is  an  integer  greater  than 
zero  but  less  than  or  equal  to  n^.  The  first  r^ 
failures  times  of  the  random  variable  are 
observed  and  recorded  as  the  failures  occur.  If 
the  first  r^  failure  times  of  the  random  variable 
Ti  are  given  by  0  <  tii  ^  t^  -  ♦  •  •  ° *■%  , 

then  the  likelihood  function  for  the  item  of  the 
ith  type  is 


a*. 


X,!  =  nil/(nt  -  rt)!  JT 

V 


[l-p(tlrl)] 


ni-  ri 


for  i-1,2,  .  .  .  ,M.  Proceeding  in  exactly  the 
same  manner  as  with  a  single  item,  the  M.L.E*  of  the 

'  9 

parameters,  ai  and  b^,  of  the  item  of  the  ith  type  are 


•  «i  =  I  $.  tijbi  +  (nl-ri)tl»1 


A 

bi 


and 


A 

bi 


•  bi  riai/  t  tij‘/J'  ln  + 


A  /l 

^ni"ri^  t ir^1  “Ln‘ tiri“^1  ln  t 


ij 
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As  before,  equations  (2$)  and  (26)  can  be  solved  by  an 

A  A 

iteration  process  for  a.±  and  bA,  for  i  =  1,  2,  .  .  .  M. 

If  all  of  the  shape  parameters  are  assumed  to  be  known, 
the  M.L.E*  of  the  parameter  a  is 

(27).  +  *lr  bl  J  /ri 

I  1 

for  i  =  1,  2,  .  .  .  ,M. 

SAMPLING  PLAN  III  -  Fixed  sample  size  with  time  truncation. 
Single  item  case  -  A  random  sample  of  size  n  is  selected 
and  each  member  of  the  sample  is  tested  for  a  specified 
time,  say  tQ*  The  time  tQ,  where  tQ  is  greater  than 
zero  but  less  than  infinity,  is  defined  as  the  truncation 
time*  The  number  of  failures  denoted  by  R  that  occur  in 
the  time  interval,  (0,to)Jare  observed  and  the  R  failure 
times  of  the  random  variable  T  recorded*  Using  this 
sampling  procedure,  R  is  of  course  a  random  variable 
which  may  take  the  value  r,  where  r  =  0,1,2,  .  .  .  n. 

If  the  observed  failure  times  are  t-^tg,  ,  *  .  ,tr, 
then  the  likelihood  function  is 

L  =  (aVv  f  (tj)  j^l-F(t0)J  nr  for  <  tQ. 

The  likelihood  function  for  this  sampling  plan  is 
formally  derived  in  Appendix  A.  Again,  substituting 


F(t)  and  f(t)  of  equations  (3)  and  (Ij.)  in  the  above 
relation,  the  M.L.E.  of  a  and  b  are  derived  in  the 
usual  manner  and  determined  to  be  given  by 


respectively.  As  before,  an  iteration  process  can  be 

A.  A 

employed  to  solve  equations  (28)  and  (^9-)  for  a  and  b. 
Also,  under  this  sampling  plan  the  M.L.E.  of  the  scale 
parameter  is  given  by 


when  the  shape  parameter  is  assumed  to  be  known. 
Complex  item  case  -  If  the  complex  item  is  assumed  to 
be  composed  of  M  individual  items,  the  experimental 
procedure  is  to  select  a  random  sample  of  size  n^  from 
the  item  of  the  ith  type,  for  i  =  1,  2,  •  .  •  ,M.  The 
sample  of  the  ith  type  is  tested  for  a  specified 
truncation  time  t*  yand  the  number  of  failures  denoted 
by  that  occur  in  the  time  interval^ (O,t^o)aare  ob- 


served  and  the  failure  times  of  the  random  variable 
recorded.  As  before,  the  number  of  failures,  R^, 
is  a  random  variable  which  may  take  the  value  r^ 

i 

where  r^  s  1,  2,  •  •  ,  ,  n ^ »  If  t^g,  *  *  #,^lr^ 

are  the  observed  values  of  the  r^  failure  times;  then 
using  the  likelihood  function. 


Li  =  ( "t)  jf*  f<tu>  [i-p(tlo)]  ni'ri  ’} 

the  M,L.E,  of  a^  and  b^^  are  determined  to  be  given  by 

•  *i  *  L  tijbi  +  ^ni“ri)ti0bi  I  /ri 


t3-i) 


A  A*  * 

(32).  =  r]£i/  L  jf,  tijbi  ln  feij 


A 

+  (n,-r,)t.  bi  In  t.  •  a.  2*  In  t 

i  i  io  xo  1  <ai 


A  A* 

Lo  *  4 

\ 

for  the  item  of  the  ith  type,  where  1=1,  2,  .  .  ,  M. 
If  all  the  shape  parameters  are  assumed  to  be  known, 
the  M.L.E,  of  the  parameter  is 


ij 


i 
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Thus,  it  is  readily  seen  that  considerable  cal- 

A  A 

culations  are  required  to  solve  for  a  and  b  in  each 
of  the  above  sampling  plans.  But  by  utilizing  a  high 
speed  electronic  computer,  the  iteration  process  re¬ 
quired  to  solve  for  a  and  ^  is  readily  accorapolished. 
Kao  £9!  outlines  a  graphical  method  of  obtaining  a 

A  A 

first  approximation  for  a  and  b.  If  these  first  ap- 

A  A 

proximations  of  a  and  b  are  provided  to  the  computer 
as  the  initial  trial,  the  iteration  method  should 

A  A 

quickly  converge  to  a  solution  for  a  and  b.  In  the 
next  chapter,  the  above  sampling  plans,  assuming  the 
parameters,  a  and  b,  are  unknown,  will  be  utilized  in 
the  proposed  models  for  analyzing  Polaris  Missile 
System  Trouble  Failure  Reports.  It  will  be  further 
assumed  that  high  speed  electronic  computers  are 
available  to  acoompolish  the  calculations  required 

A  A 

to  solve  for  a  and  b. 
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CHAPTER  IV 


APPLICATION  TO  POLARIS  MISSILE  SYSTEM 

In  this  chapter,  three  models  are  proposed  which 
utilize  Polaris  Missile  System  (PMS)  Trouble  Failure 
Report  (TFR)  as  experimental  data  for  determining 
M.L.E.  of  parameters.  By  following  the  procedures 
proposed  in  the  models  below.  It  may  be  possible 
to  obtain  a  useful  point  estimate  of  quality  indices, 
such  as  mean  time  to  failure  and  reliability,  of 
certain  items  composing  the  missile  system.  The 
models  are  applicable  primarily  to  missile  items 
for  which  a  life  history  is  maintained.  An  item 
will  be  considered  to  have  a  life  history  if: 

(1) .  It  is  possible  to  distinguish  each  item 
by  some  identification  system  from  all  other  Items 
of  the  same  type. 

(2) .  >  A  permanent  record  is  maintained  of  the 
total  operating  time  to  failure. 

(3) .  Assuming' that  the  item  doe3  not  fail,  it 

i3  possible  to  estimate  accurately  the  total  operating 
time  accumulated  by  the  item  at  some  future  date. 

In  each  of  the  models,  the  total  operating  time 
to  failure  of  an  item  is  assumed  to  be  a  random  variable 
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T  that  is  speoified  by  the  two  parameter  Weibull  prob¬ 
ability  distribution  with  both  parameters  considered 
to  be  unknown.  Accordingly,  if  n  items  of  the  same 
type  are  selected  at  random,  this  sample  is  assumed 
to  consist  of  n  independent  identically  distributed 
random  variables  whioh  are  characterized  by  the  same 
probability  law  as  T.  It  is  further  assumed  that 
none  of  the  items  have  accumulated  operating  time 
prior  to  their  purchase  by  the  United  States  Navy. 
Also,  an  item,  which  has  failed  and  been  repaired, 
is  considered  to  be  the  same  as  a  newly  purchased 
item  with  no  accumulated  operating  time.  In  each 
of  the  models,  the  operating  time  of  the  items  is 
accumulated  by  check  outs  or  other  reasons  for  op¬ 
erating  the  item,  and  the  TFR  which  provides  the 
total  operating  time  to  failure  is  used  as  the  ob¬ 
served  value  of  the  random  variable.  Naturally,  it 
is  assumed'  that  the  time  to  failure  reported  in  a 
TFR  is  accurate. 

The  first  model  proposed  is  Sampling  Plan  .1 
which  is  discussed  in  Chapter  III.  Consequently, 
for  a  sample  of  size  n,  the  M.L.E.  of  a  and  b  are 
given  by  equations  (16)  and  (17) »  respectively, 
and  a  point  estimate  for  the  reliability  and  mean 
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time  to  failure  can  be  obtained  by  substituting  the 

A  A 

calculated  values  of  a  and  b  into  equations  (5)  and 
(7),  respectively*  In  the  case  of  a  complex  item, 
the  procedure  is  the  same,  except  that  the  values 
are  obtained  for  the  M.L.E.  of  and  b^  for 

i 3 1,  2,  *  *  *  ,  M,  with  equations  (19)  and  (20)* 
Again,  a  point  estimate  of  the  reliability  and  mean 
time  to  failure  can  be  obtained  by  substituting  the 

/»  A 

values  of  a^  and  b^  into  equations  (11)  and  ( II4. ) , 
respectively*  Notice  that  equation  (lij.)  must  be 
evaluated  by  numerical  methods  to  obtain  a  point 

9 

estimate  of  X  *  The  limitation  of  this  model  is 
the  waiting  time  required  to  observe  the  nth  failure* 
Prior  to  outlining  the  next  two  models  which 
are  applications  of  Sampling  Plans  II  and  III,  it  is  . 
necessary  to  present  additional  background  concerning 
the  PMS.  Consider  an  item  of  a  given  type,  say  X, 
which  continues  to  accumulate  operating  time  by 
check  outs  until  it  fails*  Because  the  PMS  is  an 
expanding  program,  items  of  type  X  are  constantly 
being  added  to  the  system  and  also  accumulating 
operating  time*  Consequently,  if  n  items  of  type 
X  are  selected  as  a  sample,  each  of  the  items  could 
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have  accumulated  varying  amounts  of  operating  time 
depending  upon  when  the  item  was  purchased  by  the 
Navy.  The  effect  of  the  above  situation  is  to  com¬ 
plicate  item  and  time  truncation  sampling  plans  as 
discussed  in  Chapter  III.  But  Sampling  Plans  II  and 
III  can  still  be  used  to  obtain  M.L.E.  of  the  para¬ 
meters  ,if  a  life  history  is  kept  for  items  of  type 
X,  as  follows. 

The  second  method  for  obtaining  point  estimates 
of  quality  indices  employs  Sampling  Plan  II.  The 
procedure  is  to  select  P  samples  of  size  nj,  where 
Jsl,  2,  .  .  •  ,P,  from  items  of  type  X.  The  items 
in  the  jth  sample  are  selected  so  that  each  of  the 
members  of  this  sample  are  expected  to  have  accumu¬ 
lated  approximately  the  same  amount  of  operating 
time  at  a  future  date.  For  each  of  the  P  samples, 
fix  in  advance  the  number  of  observed  failures,  de¬ 
noted  by  rj,  which  will  terminate  the  jth  experiment. 
If  equation  (21)  is  used  as  the  likelihood  function 
for  the  jth  sample  of  3ize  n^  and  If  the  P  samples 
are  assumed  to  be  independent,  the  likelihood 
function,  denoted  by  £.  ,  of  the  P  samples  is  given  by 

£  rz  rfrn}l/(nrT})  TT  fttp 
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Upon  substituting  F(t)  and  f(t)  of  equations  (3)  and 
(4)  in  the  above  relation,  the  M.L.E.  of  the  para¬ 
meters,  a  and  b,  are  determined  in  the  usual  manner 
to  be  given  by 


!' 


Therefore,  a  point  estimate  of  a  quality  index  can  be 
obtained  by  substituting  equations  (34)  and  (35)  into 
the  equation  which  defines  the  quality  index  desired. 

In  the  case  of  a  complex  item,  the  method  is  completely 
analogous  except  there  are  M  items  of  different  types. 
Again,  life  histories  must  be  maintained  for  each  type 
of  item.  > 

The  third  model  uses  Sampling  Plan  III.  As  in  • 
the  second  model,  P  samples  of  size  n^,  where 
J  -1,  2,  •  •  •  ,P,  are  selected  from  items  of  type 
X  so  that  each  member  of  the  jth  sample  is  expected 
to  have  accumulated  the  same  amount  of  operating 
time  at  a  future  date.  A  truncation  time,  denoted 
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by  tjo»  is  selected  for  each  of  the  j  samples  and  the 

experiment  is  considered  to  be  terminated  after  the 

truncation  time  has  elasped  for  each  of  the  -j  samples 

of  size  n^,  where  j  -1,  2,  •  •  ,  ,P,  As  in  Sampling 

Plan  III,  the  number  of  failures,  denoted  by  Rj,  which 

occur  during  the  truncation  time,  t  *  ,  in  the  jth 

Jo 

sample  is  a  random  variable.  Assuming  that  each  of  the 
P  samples  are  independent,  the  likelihood  function,  de¬ 
noted  by  o£  ,  is 

•  £  -  7T  l, 

J 


where  Lj  is  the  likelihood  function  of  Sampling 
Plan  III,  Proceeding  as  in  Sampling  Plan  III, 
the  M.L.E.  of  a  and  b  are  given  by 


and 


+  lW* 
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In  the  case  of  a  complex  item  the  procedure  is  com¬ 
pletely  analogous.  Therefore,  point  estimates  of 
quality  indices  can  be  determined  in  both  the  single 
and  complex  item  case  by  substituting  the  equations 
for  the  M.L.E.  of  a  and  b  into  the  equation  for  the 
desired  quality  index. 

In  the  last  two  models,  the  selection  of  the 
members  of  the  jth  sample  should  occur  naturally 
since  a  group  of  items  of  a  given  type  is  usually 
issued  to  a  submarine  at  the  same  time.  These 
items  should  accumulate  operating  time  at  the  same 
rate  through  cheok  outs  until  failure.  The  time 
when  each  of  the  items  is  issued  can  be  determined 
by  examining  the  life  history  of  the  items;  con¬ 
sequently,  the  members  of  the  jth  sample  can  be 
determined  for  either  model  two  or  three. 
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APPENDIX  A 


In  Sampling  Plan  III,  n  items  are  tested  until  a 
specified  time,  tQ,  has  elasped.  The  number  of  failures 
denoted  by  R  that  occur  in  the  time  interval  (0,to)  are 
observed  and  the  R  failure  times  recorded.  Using  this 
sampling  procedure,  R  is  a  random  variable  which  may 
take  the  value,  r,  where  r  =  0,  1,  •  .  •  ,  n.  As 
before,  equations  (3)  and  (Ij.)  define  the  probability 
density  and  cumulative  distribution  functions  res¬ 
pectively  of  the  random  variable  T.  The  discussion 
which  follows  i 3  a  formal  derivative  of  the  likelihood 
function  for  Sampling  Plan  III. 

Maximum  likelihood  estimators  are  functions  of 
the  data  obtained  from  an  experiment.  In  this  case, 
the  experiment  is  to  specify  a  truncation  time  tQ 
and  observe  the  number  of  failures,  R,  and  the 
failure  times,  T^,  which  occur  in  the  time  in¬ 
terval  (0,to);  consequently,  the  joint  density 
function  of  the  number  of  failures,  R,  and  the 
times  to  failure,  T^,  given  that  the  failures 
occur  in  the  time  interval  (0,to)  would  be  ap¬ 
propriate  for  determining  the  likelihood  function. 
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Let  A  be  defined  as  the  event, 

T2<  ^2»  •  •  •  *  TR<  tjj)  where  the 
are  not  ordered  and  t±<  tQ  for  i  = 1,  2,  .  .  ,,r. 

Let  B  be  defined  as  the  event,  (R=r),  where 

***  0,  1,  .  .  .  ,  n. 

Let  C  be  defined  as  the  event  that  the  failures 
which  occur  do  so  before  time  t0. 

Let  the  event,  AB,  be  expressed  as  (T1<  tQ, . . . ,Tr < tQ) 
Then  the  problem  is  to  obtain  the  distribution  function 
which  represents  the  probability  statement  P(AB|c). 

Using  conditional  probability  relations,  it  is  ob¬ 
served  that 

P(ABlC)  -  P  (  ABC )  ,  P(BC)  _  P(AlBC)P(B|C) 
pTbcT  P(C)  " 

where 

P(AlBC)=  P(T1<t1#  •  .  .  ,TR<tR|Bnc) 

tl*  *  *  *  »Tr<tr,Tr.j.j>to> * « « ^n^o  B^J 

P(BHC) 
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and 


P(Bjc)  =  gJiUoQ'  jJ-P(t0)j  n"r 

therefore 

P(ti,t2M..,tr,Rsr)s  jJ-F(t0)j  n"r 

for  t±-ft0 

Accordingly,  the  joint  density  function  is  obtained 
by  taking  the  partial  derivatives  of  the  above  dis¬ 
tribution  function  with  respect  to  t^  for  i=l,*..,r 
and  is  given  formally  by 


f (ti, • . • , tr,R s  r) 


oaj*(ti>  n'r 

for  t^<.  tQ 


38 


